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Abstract 



A simple derivation of the Jalilian-Marian, lancu, McLerran, Weigert, 
Leonidov and Kovner (JIMWLK) equation for the evolution of small-x QCD 
wavefunctions is given. The derivation makes use of the equivalence between 
the evolution of a (in general complicated) small x wavefunction with that 
of the evolution of the (simple) dipole probing the wavefunction in a high 
energy scattering. 

1 Introduction 

One expects the light-cone wavefunction of a very high energy hadron to 
exhibit saturationjQ. That is, the occupation number, in transverse phase 
space, of quarks and antiquarks of momentum fraction x should be near 1 
for transverse momentum k± < Qs{x) while that of gluons should be of size 
1/a. Qs, the saturation momentum, depends on the hadron in question and 
on the longitudinal momentum fraction of the parton being considered. Qs 
values of about IGeV, or perhaps a bit larger, may be available for study at 
HERA and at RHIC and it is extremely pressing to find out if saturation is 
clearly visible at these accelerators. 

Over the past 7 years or so there has been an ambitious program dedicated 
to finding appropriate equations for dealing with such high density QCD 
wavefunctions[|, ||, |, ||, 0. This program has been quite successful and 
a renormalization group equation in the form of a functional Fokker-Planck 
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equation for the wavefunction of a high energy hadron has been given by the 
authors of Refs.(4-6) (JIMWLK). The most complete derivation is given in 
Ref.(6) where the equation is written in terms of a covariant gauge poten- 
tial, a, coming from light-cone gauge quanta in a high energy hadron. The 
mixture of light-cone gauge and covariant gauge comes about as a mathemat- 
ical technique for simplifying an equation which appears to be exceedingly 
complicated when expresssed purely in terms of light-cone gauge quantities. 

It has been known for sometime that there is a dual description where 
the unitarity limit for scattering of a small dipole on a high energy hadron is 
equivalent to saturation[|, |[. In this dual picture one views dipole-hadron 
scattering in a frame, the dipole frame, where the dipole is left-moving and 
the hadron is right-moving. The dipole rapidity, y, is kept small enough, 
« 1, so that one need not consider radiative corrections to the dipole 
wavefunction. In a high energy scattering the rapidity of the hadron is very 
large and the hadrons wavefunction involves many gluons. If, at a given 
impact parameter, the S-matrix for dipole-hadron scattering is significantly 
less than 1 then this strong scattering is equivalent to the statment that 
the occupation number of quarks and antiquarks, at a transverse momentum 
corrresponding to the inverse dipole size, is of the order of 1. 

In this paper we use the dipole frame to derive the JIMWLK equation 
in a rather simple manner. As one increases the energy of a high energy 
dipole-hadron scattering by a small amount, this increase can either be put 
into the hadron or into the dipole. If the increase is put into the hadron 
then the wavefunction of the hadron changes (evolves) . If the increase is put 
into the dipole the dipole can emit a gluon in an easily calculable manner. 
By equating these two pictures one arrives at the JIMWLK equation. In a 
sense our physical picture is a simple way of viewing the set of mathematical 
transformations carried out in Ref.(6) in their derivation of this equation. 

2 The JIMWLK Equation 

Let 

/oo 
dx^A+{x,xJ)} (1) 
-co 

be a path-ordered line integral taken, say, at x+ = where 
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= Y: T^Alix) (2) 

a=l 

with T'* the generators of color SU{Nc). The symbol P in (1) indicates that 
matrices A+{x,X-) having larger values of x_ are put to the left of those 
having smaller values of X-. In a moment we shall specialize our discussion 
to Coulomb gauge, but for now any gauge for which is not identically zero 
and for which Ai{x, X-) — * is suitable. For simplicity of presentation we 

take the T" to be in the fundamental representation of SU{Nc). 

Let 0{A^) be a gauge invariant functional formed out of V's and V^^'s at 
various values of x. For example, 

02{x,,X2)=trV\x,)V{x2) (3) 

and 

04(xi,X2,X3,X4) =trV\xi)V{x2)V\x3)V{x^) (4) 

constitute such functionals so long as Ai vanishes at X- = ±oo. 
Suppose \p > is the state of a hadron having momentum 

{po,p,Pz) = M{cosh Y, 0, sink Y) (5) 
where the rapidity Y is large and positive. Then one may write]^] 

< O >Y=< p\0\p >= J D[a{x,x_)]0{a)WY[a] (6) 
where the weight function Wy is given by 

WY[a] = I D[A^,]6{A+ - a)6{F{A))Ap[A]e'''^^\ (7) 

The symbol D indicates a functional integral, F is a gauge fixing, and Ap 
is the corresponding Fadeev-Popov determinant times an operator which 
projects out the state \p > . We suppose Wy, the distribution of values 
of the potential A^ in the hadron wavefunction, is normalized to 

J D[a]WY[a] = 1 (8) 
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To interpret (6) more completely we now specify our gauge choices as 
Coulomb gauge. The advantage of Coulomb gauge[|l^] is that it gives a nice 



separation between left-moving quanta, right-moving quanta and gluons ex- 
changed between left and right-movers. Quanta having k^/k_ » appear 
as they would in = gauge. Quanta having k^/k^ >> 1 appear as they 
would in A_ = gauge. Exchanged gluons appear as in covariant gauge if 

/c_ << k±. Thus < 0{xi,X2) >Y corresponds to the S-matrix for scatter- 
ing of a left-moving dipole, of size Xi — X2 on a right-moving hadron. Wy [a] 
gives the field = a coming from all the right-moving quanta in the high 
energy hadron. The rapidity of the scattering is Y. This latter statement is, 
perhaps, not immediately clear since O2 is defined in terms of left-movers ex- 
actly on the light cone x+ = 0. However, in order that O2 correctly represent 
a left-moving dipole we must assume that its rapidity, y, obeys agi/ « 1. 
If c^sH > 1 the bare dipole will become dressed with additional left-moving 
gluons so that O2 would not represent a physical dipole. 

Now consider the y-dependence of < O >y . From (6) one can write 



^<0>y= J D[a]0{a)^Wy[a]. (9) 



However, one can equally well imagine calculating ^ < O >y by evaluating 
the change of the left-moving dipole, due to an additional gluon emission, as 
one increases the rapidity of the left-moving system by an amount dY. The 
change in the left-moving state is given by the gluon emission and absorptions 
shown in Figs.l and 2 where spectator quark and antiquark lines are not 
shown. In Fig.l we have assumed that a gluon line connects a quark and an 
antiquark. This is for the sake of definiteness. We could have equally well 
assumed a connection between two quarks or between two antiquarks. In 
Fig. 2 we suppose that the gluon is emitted and absorbed on the same quark 
(or antiquark) line. The vertical lines in the figures represent the "time," 
a;_ = 0, at which the left-moving system passes over the right-moving hadron 
and its gluon fields. This view of the F— dependence of < O >y in terms of 
a change of the rapidity of the (rather simple) left-movers is in the spirit of 
work previously done by Balitsky|Tl[], Kovchegov|jT2[], and Weigert||H]. 



Let's examine the contribution of the graph of Fig. la to ^ < O >y in 
more detail. In this contribution a gluon is emitted off a quark at x^ and 
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x_ 

(a) (b) (c) (d) 

FIG.l 

absorbed by an antiquark at The emission and absorption happen after 
the left-moving system has passed the hadron so the gluon does not interact 
with the hadron. We do the calculation in = light-cone gauge which 
for left-movers is equivalent to Coulomb gauge. Including the term where 
the gluon is first emitted off the antiquark and then absorbed by the quark 
one has the simple replacement [|13] . 



in going from O to -^O. The notation ®V means that the color indices 
on V and V"^ are fixed with no matrix multiplication implied. Of course 
T'^V^{xi) does involve a matrix multipication of T'^ times V^^Xi) on the 
right-hand side of (10). The fact that comes to the left of V^Xi) but to 
the right of V{x2) (10) follows from the a;_— ordering given in (1) and 
the fact that the contribution of Fig. la has gluon emission and absorption 
at (large) positive values of x_, that is after the scattering has taken place. 

Now consider the contribution illustrated in Fig.lc. This contribution 
involves the replacement 



TT- J [Xi - Z)^[X2 - Z) 

(11) 

The new elements which arise in (11) compared to (10) are the following: 
(i) The overall sign is (-) because the gluon has been emitted before the 
scattering and absorbed after the scattering, (ii) The factor Vcd{z) gives the 
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interaction of the gluon with the right-moving hadron as an eikonal "phase" 
as in (1) but now in the adjoint representation, (iii) Now the T*^ come to 
the right of V^(xi) because the gluon is emitted off the Xi— line before the 
interaction with the hadron. For reason which should become clear shortly 
it is useful to write 

V\x,)T'^ = V\x,)T'V{x,)V\x,) = VUx,)T''V\x,) (12) 

in order to bring the T to the left of V"^(xi). Thus one can write, replacing 
(11), 

a /■ ,2 (^1 - ^) ■ (^2 - z) 



■ [VKx^V{z)UT''V\x^ ® V{x,)T'. (13) 

There is a simple mnemonic for seeing the structure of the V's and V's in 
(13). Instead of the graph in Fig.lc we could use that shown in Fig.3 where 
an adjoint line integral starts at large positive values of x_ follows the quark 
line at decreases to large negative values and then proceeds from 

large negative to large positive values of x_ at the point z_. The interaction 
of the hadron with the quark and gluon, both at Xi, in Fig.3 is the same as 
with the quark line in Fig.lc. 

It is now straightforward to add in the contributions of the graphs shown 
in Fig. lb and Id to get the total contribution of the graphs of Fig.l as 



vr^ J (jxi — z)'^{x2 — z)'^ 



■{l + V\x,)V{x2)-V\x,)Viz)-V\z)V{x2n (14) 
when O ^ 47O. This can be written as 

dY 



withfl 
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X 



(a) 



(b) 
FIG.2 



X 

(c) 



^ J {2t:Y {x- zY{y- zf 

■ {1 + V\x)V{y) - V\x)V{z) - V\z)Vy)}ab (16) 

where it is understood that the two functional derivatives on the right-hand 
side of (15) do not both act on cither V or on 

One might expect that the contributions of the graphs shown in Fig. 2 
would come from an equation like (15) but where the two functional deriva- 
tives act on a single term, say V\x]), and where the other factor, ^(^2) is 
absent. This is almost the case. The contributions illustrated in Figs. 2a and 
2b do correspond to an equation like (15) taking only the first two terms in 
the brackets for the r] given in (16). However, 

■ V\b) = ^ / j£^,{V\x,)vHz) + VHz)V{x,)UT'T'vHx,) (17) 
while the contribution of the term shown in Fig. 2c is 



(18) 
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FIG.3 



The difference between the expressions on the right-hand sides of (18) and 
(17) is 

= - — y j-—-^^tr[rV\x,)V{zWV\x^ (19) 
Thus taking (15) and (19) we can write 



(20) 

where 

9< = 2^/ ^£^,tr[T^V\x,)V{z)] (21) 

and where, as usual, x_ and y_ are positive. 

Finally, equating these two pictures of F— evolution given by (9) and by 
(20) we can write 

j D[a]0{a)^[a] = J D[a]0{a) 

(22) 



where we have integrated by parts in going from (20) to the right-hand side 
of (22). 

Eq.(22) can be cast into a form closer to that used in JIMWLK by defin- 
ing, for F — F >> 1, 

Wyyia] = Wvia] (23) 

where 

a{x,x-) = e^'^aix.e^'^^x-). (24) 

That is, VFf> y[a] corresponds to the same functional integral given in (7) but 
where the state \p > now has rapidity Y and where only fields of longitudinal 
momentum k^/p^ > are included in the functional ingtegration. Since 
the right-moving fields are effectively in = light-cone gauge in (7) the 
Lorentz boost indicated in (23) and (24) is purely kinematic. Since 0(a) in 
(22) constitutes a very general form of gauge invariant functional we presume, 
but do not know how to prove, that (22) - (24) imply 

the form of the JIMWLK equation given in Refs. 6 and 7. 

What are the values of x_ and y- in (25), or equivalently, in (22)? The 
derivation of (22) makes clear that all that is required is that x_ and y- be 
large enough so that 

g / a{x,x'_)dx'_, g / a{y,y'__)dy'^ « 1 (26) 

allowing (14) to be replaced by (15). The exact values of x^ and ?/_ are not 
important although one should require]^, |^ 

^e-(^-^) « i/_ « ^e-^y-y-^y) (27) 

with Ay = 1/a in order that there not be additional gluons emitted with 
rapidity between Y — Y and the rapidity of the gluon we are considering. 
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Q is a scale characterizing the typical transverse momentum of gluons at 
rapidity Y — In terms of kinematics of the original frame used in (22) 
the requirement (27) is simply that as(in{x^Q) « 1 and asin{y^Q) « 1 
similar to the condition imposed earlier in the rapidity of the left-moving 
dipole. 

We have noted several times that the right-moving gluons determining 
W can be viewed as light-cone gauge, = 0, quanta. Operations, such as 
indicated in (12), which bring generators always to the left of and to the 
right of V make the further choice of an ^4+ = gauge where the gauge fields 
vanish for x_ < 0, thus matching the choice of Refs.6 and 7. Equally well one 
could have chosen to move all generators, T°, to the right of V"^ and to the 
left of V which would correspond to the light-cone gauge where potentials 
vanish, for high momentum particles, for x_ > as was done previously in 
Refs. 14-16 in a different context. 
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